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Abstract 

The jets are the final state manifestation of the hard parton scattering. Since at LHC 
energies the production of hard processes in proton-proton collisions will be copious and varied, 
it is important to develop methods to identify them through the study of their final states. In 
the present work we describe a method based on the use of some shape variables to discriminate 
events according their topologies. A very attractive feature of this analysis is the possibility 
of using the tracking information of the TPC+ITS in order to identify specific events like jets. 
Through the correlation between the quantities: thrust and recoil, calculated in minimum bias 
simulations of proton-proton collisions at 10 TeV, we show the sensitivity of the method to 
select specific topologies and high multiplicity. The presented results were obtained both at 
level generator and after reconstruction. It remains that with any kind of jet reconstruction 
algorithm one will confronted in general with overlapping jets. The present method determines 
areas where one does encounter special topologies of jets in an event. The aim is not to supplant 
the usual jet reconstruction algorithms, but rather to allow an easy selection of events allowing 
then the application of algorithms. 
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Figure 1: Plot: 77 vs. A(<f> — 7r/2) vs. pt for a typical event. The leading particle is placed at ir/2 
rad. Different radius are shown. 



1 Introduction 

In QCD the jets are defined as cascades of consecutive emissions of partons initiated by partons 
from any initial hard process. The partons produce the observed hadrons due to confinement. 
Di-jets were discovered in 1975 in e + e~ collisions [I], the observation of three coplanar jets 
has provided the first experimental evidence for the existence of gluon [2] |3j [4] [5] . 

Event shapes measure geometrical properties of the energy flow in QCD states. Especially 
in the case of e + e~ collisions and Deep Inelastic Scattering, they were among the most studied 
QCD observables, both theoretically and experimentally. Sphericity was used at SLAC to show 
the evidence for the existence of jets in the annihilation process: e + e~ at energies up to 7.4 
GeV in the c. m. [T|. In 1979 the collaboration MARK- J used the variable oblateness for 
describing processes where three prolong jets were produced at energies up to 31.5 GeV in the 
c. m. [2]. 

Of course the shape variables are calculated in terms of the final partons and they are 
measured in terms of the final observed hadrons, so it is possible to use them in order to do 
studies of corrections by hadronization effects [S] . Also a vast number of strong coupling 
constant measurements was reported [6|- In the case of the MC generators the validation of 
some results were also published [101 111] . 

Experimentally, jets in proton-proton collisions are defined as an excess of transverse energy 
over the background of the underlying event^l with a typical cone radius Rc ~ 1 in r\ — <f>. 
In ALICE we do not yet have an extensive calorimeter. For this reason, in order to define 
and reconstruct jets we are using tracking measurements. The jets identification is based on 
different algorithms; for example the cone ones. In this work we propose a method which uses 
the event shapes to identify interesting topologies of events like jets, which of course does not 
prevent the later use of jet algorithms on the selected events. 

The main problem we will face at LHC is the big number of jets overlapping one another. 
This makes all analysis subject to numerous cuts which do increase the systematic uncertainties 
of the results. A typical events represented in the eta-phi plane is shown in Fig. [T] To guide 
the eye the contours of jets with _R cone =0.3, 0.5 and 1 are shown for the ALICE acceptance. 
However in this note we would like to emphasize the fact that there is a possibility, using the 
event shapes, to identify painlessly events where the topology is much simpler than shown in 
Fig. [T] In the present note we first describe the technicalities of the method of event shape 
analysis giving some overview of salient topologies like events where the jets are only partially 
detected in the acceptance, events with 2 jets and events with 3 jets. Then, the main features 
of each class are studied and demonstrated, like total pt spectra of jets, the sensitivity on 
the parameters variation within the method, the yield of specific topologies in function of 
multiplicity. Finally we introduce also a bulk analysis of all events without selection. We used 
this to present a comparison between Pythia and PHOJET events for different parameter 
values. 



1 Underlying events are formed from the beam-beam remnants, initial-state radiation and possible from soft and 
semi-hard multiple parton interactions. 




Figure 2: Shape variables at hadron colliders are defined over particles within a central region, C, 
and by the addition of a recoil term also defined within C, are sensitive to the forward region C. 



2 Definitions of the shape variables 

By construction the following quantities are Infra Red and Collinear (IR&C) safe. The shape 
variables at hadron colliders are defined over particles within the acceptance of the detector. 
The thrust (T) is defined as in the e + e~ case, but using only transverse variables [12] : 

rsmM %^ (1) 

n t 

In the literature it is more common to find the following definition: 

r= 1 — T (2) 

(related to the sphericity of the event). 

The range of r is between (for events with narrow back-to-back jets) and 1/2 (for events 
with a uniform distribution of momentum). 

The recoil term 7? is simply the vector sum of the transverse momentum: (Fig. [2} 

^cki 1 ? 7 " 1 (3) 

This quantity measures the balance of momenta of the event. For example for a di-jet event, 
with only one jet inside the acceptance of the detector (monojet in the further text) : R tends to 
1, because there are no vectorial cancellations in the numerator which appear in the definition 
of R. Otherwise, in the case of the perfect 2 back-to-back jet completely inside the acceptance: 
R tends to 0. 



3 Description of the event shape analysis (ESA) 

We have computed the shape variables at level generator and reconstruction. The partici- 
pants in the computations are primary charged particles in the generation case, and tracks 
associated with primary particles in the reconstruction one. In both cases events are selected 
using the trigger criteria and we restrict the analysis for events with primary vertex in 

z direction: \v z \ < 10 cm. 

The requirements to perform the computation are: 

1. Generation level: 

Event level: the first step is for selecting hard events: p \ eadin 3 > 3 QeV/c and \ v leadin v\ < 0.5. 

This requirement guarantees that the leading jets are contained within the ALICE 
acceptance with a AR cone = \J Arf + A(f> 2 < 0.7 

2 Primary particles are defined as particles produced in the collision, including products of strong and electro- 
magnetic decays, but excluding feed-down products from strange weak decays and particles produced in secondary 
interactions. In the simulation these are the final state particles created by the event generator, which are then 
propagated (and decayed) by the subsequent detector simulation. 

3 This trigger uses the logical AND between the V0-OR signal and the Pixel-Fast-OR. And it is an option if one 
needs to assign each trigger to a specific bunch crossing |14| . 
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Figure 3: Illustration of the distribution of the momenta in the transverse plane for two classes of 
events, (left) Dijet event, (right) monojet event. In both cases the thrust axis is very close to the 
direction of the leading particle. 



Particle level: for primary charged particles in the acceptance: \q\ < 1 and p t > 1.5 GeV/c, the 
shape variables are computed. The cut in p t eliminates the mostly soft underlying 
component. 

2. Reconstruction level: In the present analysis we used tracks reconstructed by the 
TPC and ITS. The event cuts described above are applied. 

Track level: tracks associated to primary particles in the acceptance: |7;| < 1 and with p t > 1.5 
GeV/c. To select this class of tracks we applied the following cuts. 

(a) TPC refit. 

(b) At least 50 clusters in TPC. 

(c) Covariance matrix cuts. 

(d) Reject kink daughters. 

(e) Maximum DCA (in xy and z) to vertex 3 cm. 

For more details see [13] . 

In the Fig. [3] we show a picture of two classes of candidate events: di-jet and mono-jet (to 
separate these events we used their values of shape variables according with the discussion 
in the next paragraphs). For a given event we have taken the projection of the particles 
momentum in the transverse plane. It is well visible that the thrust axis is very close to the 
direction of the transverse momenta of the leading particle (particle with the highest p t in the 
event). 

The method starts by plotting a two dimensional distribution ("thrust map"), with r 
(1 — T) in the horizontal axis and R in the vertical axis. This plot allows to identify different 
classes of events according with their location in the thrust map. 

1. Region A. Suppose a di-jet event which occurs completely inside the ALICE acceptance 
(M < !)■ I n this case, we have in the transverse plane; the thrust axis [Tit) almost 
collinear to the direction of the leading particle. So, the ratio of definition (1) tends to 
1, and the recoil term tends to 0. So, the region A is characterized by events with small 
values of 1 — T and R, corresponding to di-jet events. 

2. Region B. Events with only one jet in the acceptance of the detector will have small 
values of r and due to the absence of vectorial cancellations in the numerator of the 
recoil term, these classes of events will have the biggest values R. So, the region B is 
populated by monojet events. 

3. Region C. The most isotropic events (with high 1 — T and small R) of the sample have 
to compensate the transverse momentum. This zone of the thrust map is characterized 
by the presence of three jet events that we like to call, as at LEP, "mercedes" events. 

The intermediate region between A and B is populated by the combination of mono-jets 
and "incomplete di-jets" ; di-jets are incomplete due to their high R, many particles are outside 
of the acceptance. 



The intermediate region between B and C is populated by events were 3 or more jets are 
emitted but at small angles with respect to the direction of the leading jet. 

In the present note we will not deal with these cases which will be the subject of a special 
study at later time. 



4 Event shape analysis in minimum bias simulations. 

The present analysis uses the standard simulated samples staged on alicecaQ. The production 
corresponds to PDC09, minimum bias simulations of proton-proton collisions at 10 TeV in the 
center of mass, the events were generated with Pythia. The simulation of the detector included 
a magnetic field of 0.5 T. The results were obtained through the analysis of 1, 200, 000 events. 
In the figure [4] the thrust map is shown for generated (upper) and reconstructed (bottom) 
data. In both histograms two interesting regions are exhibited (A and B regions described 
before). Note, that any exceptional event with an azimuthally uniform distribution would be 
immediately "detected" in the unpopulated area: 1 — T > 0.35. 

The line which appears in the high R part of the map is due to the definition of the 
variables. For example, suppose that we have one event with a certain thrust axis nto, and 
in the event there are N particles. If the vectorial transverse momentum of the particle i is 
p t,i, then from the definitions of R and T: 

l- r + fl= l- ^lf^- 1 f»l + ^4 ( 4 ) 



Now, we use the fact that: 



Then: 



\^2^pu\ < I"pmI ( 5 ) 



1-T + R<1-T+1 (6) 
The maximum value which T can reach is 1, so we find the following restriction: 

1 - T + R < 1 (7) 

In the Fig. [5]there is a plot which shows the correlation of the 1 — T values computed from 
the generator information and from the reconstruction. Note that there is a small leakage 
of events with 1 — T mc > 0.03 to the reconstructed zone associated to events with 2 back- 
to-back jet structure. This can be understood in terms of the reconstruction effects. For 
example a generated event with three primary charged particles (p t > 1.5 GeV/c) distributed 
isotropically in the azimuth could be reconstructed as a event with a structure of two jets. 

In the Fig. [6] we show the analogous plot for the variable R. 

For describing the topology of the events we performed an azimuthal correlation. The idea 
is to select the leading particle, and apply a rotation placing it at n/2 rad. After that, we 
plot the azimuthal distribution of the associated particles with respect to the leading one. In 
the following our convention about the azimuthal correlation it will be referred as Acj> and it 
refers to: A(cf> — 7r/2). 

If we concentrate ourselves to events in the region 1 — T < 0.03, and small R < 0.35 values. 
The azimuthal correlation of the Fig. [7] (right panel) shows that the width of the away side peak 
does not change if we modify the cut in the R range. This is in agrement with our assumption 
which suggests that R is important to select complete events in the acceptance. On the other 
hand, if we select events with small value of R (R < 0.35) and modify the range of 1 — T (Fig. 
[7] (left panel)) there is a clear evolution in the structure of the away side peak. If we increase 
the value of 1 — T the selected events include different configuration multijets (split jets) which 
are manifested in the azimuthal distribution behavior. This suggests that the cut 1 — T < 0.03 
is almost equivalent to select particles inside a cone radii R c = y (A<^>) 2 + (A77) 2 = 1. 

In the following we will concentrate on particular parts of the map, plotting the azimuthal 
correlations encountered. 

The procedure is: 

• Select events according to their values of shape variables as shown in table [T] 



4 The CERN Analysis Facilities for ALICE (alicecaf) is a cluster at CERN running PROOF (Parallel ROOT 
Facility) which allows interactive parallel analysis on a local cluster. 
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Figure 4: Thrust map: a) generator level (upper), in this case we use primary charged MC particles; 
and b) reconstruction level (bottom), in this case we used tracks associated to primaries from 
TPC+ITS information. 




Figure 5: Correlation between the quantity 1 — T computed from MC information and from 
TPC+ITS information. 




Figure 6: Correlation between the quantity R computed from MC information and from TPC+ITS 
information. 




Figure 7: Sensitivity of the azimuthal correlation of the away side peak on the thrust variable 
range. Fixed 1 — T interval and variation of the R range (right panel). Fixed R and modification 
of the range 1 — T (left panel) 



Region 


Kind of event 


Variables 


A 


Dijets 


R < 0.35, t < 0.03 


B 


Mono jets 


R > 0.9, r < 0.03 


C 


Mercedes 


R < 0.4, r > 0.25 



Table 1: 1-T and R parameters used for the present analysis. 
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Figure 8: Azimuthal correlation for particles with pt > 1.5 GeV/c: dijets (red-region A), monojets 
(blue-region B) and the mercedes events (green- region C). True spectrum (solid line), measured 
spectrum (dotted line). The leading particle is not shown. 



The azimuthal correlations of events sited at different regions of the thrust map are plotted 
in Fig. [8] The evolution of the away peak (formed by particles with it < Acj> < 2n) is 
interesting. 

As predicted, the events of region B really have a monojet topology in the azimuth. One 
can see there are associated particles which go near to the leading particle (peak in the toward 
side: < A<f> < n) but in the away side there is no corresponding jet. In contrast, for events of 
the region A the peak of the away side is located at A(j> ~ Sn/2, so, in the transverse plane we 
have 2 back-to-back jets. In the case of the events of region C we found a three-jet structure, 
in the green distribution we observe three peaks in the spectrum: the first (associated to 
the leading jet) at A<f> ~ tt/2 and the others at A(j> ~ 7n/6 and A<j> ~ ll7r/6 respectively. 
Due to the observed topology we named the latter "mercedes" events. The result bears some 
resemblance with the away side structure observed at RHIC e.g. PHENIX collaboration in 
heavy-ion collisions [15] . This observation brought us to study the presence of the same double 
hump structure at RHIC energies [TS] 

In the Fig. [9] we show for the different regions of the thrust map, a two dimensional 
distribution: pt vs. A(f> for the associated particles. 

The table [5] shows a summary of the analyzed events. 

The efficiency represents the number of reconstructed events of given topology with respect 
to the generated events. 

One observes a ~ 70% efficiency with respect to the MC. This may indicate that due to 
reconstruction and/or absorption of particles one relocates events in other parts of the thrust 
map as shown in Figs. 15161 
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Figure 9: Transverse momentum vs. A(<j> — ir/2) for the associated particles: dijets (top), monojets 
(middle) and the mercedes events (bottom). True distribution (left plots), measured distribution 
(right plots). The leading particle is shown, and you can note the cut pt > 1.5 GeV/c which we 
imposed. 



Event 


MCtracks 


ESDtracks 


Efficiency 


1 — T (t) cuts 


R cuts 


With T 


28920 


24960 


86% 


no 


no 


Dijet 


1503 


1316 


67% 


r < 0.03 


R < 0.35 


Mono jet 


8903 


8329 


79% 


r < 0.03 


R > 0.9 


Mercedes 


523 


439 


62% 


r > 0.25 


R < 0.4 



Table 2: Number of events with particular topologies in the analyzed sample of 1.2 million of 
events. 




Figure 10: Visualization of one dijet event from the region A. The lines are primary monte carlo 
tracks with p t > 0.5 GeV/c. The superimposed arrows separate each jet for a better visualization. 



Note that, out of the 1200000 events only 2.4% (with given T) of them pass the cuts imposed 
(at least 1 particle with p t > 3 GeV/c and \ r f eadln 9\ < q 5). -phe dijets in acceptance reach 
about ~ 0.125% of the total, and about 0.04% belong to mercedes event types. About 37.8% 
of the accepted events belong to clearly identifiable categories, while the others correspond to 
events with multi-jets closer to the leading jet. 

In the next section we will study each one of the topologies selected. That study includes 
a visualization of the events through the use of the tools of AliRoot[3 

4.1 Dijets and monojets 

In the present section we are investigating the p t spectra and the multiplicities of the jets in 
various configurations. As a first step we turn to the visualization tool of events in ALICE. 
We selected events located in the region A and scanned them. Due to their small R values, 
they should be inside the acceptance. Fig. 1101 shows the visualization of one of them. This 
particular event has R — 0.00132 and 1 — T = 0.0000896. The lines are primary monte carlo 
tracks with pt > 0.5 GeV/c. The figure clearly shows the whole event contained inside the 
TPC. 

Looking at more events shows always the same structure in the visualization. 

Further we computed the total transverse momenta in each jet (here we used all events of 
the region A). In order to do it, we divided the event into two parts: toward and away. The 
first (near side) contains the leading particle and it is formed by all particles in the interval: 
7r/4 rad < A(f> < 37r/4 rad. The away side is formed by particles in the interval: 57r/4 rad 
< A<j}< 7tt/4. 

The Fig. II II shows the distribution of the transverse momentum spectrum of each jet (sum 
in each azimuthal region of the transverse momentum of all the participants with pt > 0.3 
GeV/c) for generation and reconstruction. The away side distribution is left shifted ~ 1 
GeV/c with respect to the near side one. This can be understood in terms of the fluctuation 
of the neutral component of the associated jet and also as effects from the acceptance in the 

5 AliRoot is the name of ALICE Off-line framework for simulation, reconstruction and analysis. It uses the ROOT 
system as a foundation on which the framework and all applications are built. 
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Figure 11: Distribution of the sum of the transverse momenta of particles with p t > 0.3 GeV/c 
in the toward and away regions for dijet events. The toward region (black, leading jet) is formed 
by primary charged particles in the azimuthal range: 7r/4 rad < A(f> < 37r/4 rad while the away 
region (red line) corresponds to particles with: 57r/4 rad < A(j> < 7ir/4 rad. The plot shows two 
cases: reconstruction (dotted line) and generation (solid line). 



associated jet. In order to illustrate this arguments we plotted the distribution of the ratio: 
total transverse momenta in the "toward" side over the total transverse momenta in the "away" 
side for the generation case. The Fig. [12] (right panel) shows the behavior of such distribution 
(red line). The distribution manifest a clear peak at 1, this fact is in agrement with our 
assumption about the dijet structure. However the distribution shows many events where the 
near side jet represent up to 5 times the energy of the away side jet. In this respect the role of 
the value R is important, because as you can see the width of the distribution decreases as the 
R interval is decreased. The same analysis can be performed including the neutral component 
(left panel) . The away side jets with a low transverse momentum correspond to events where 
the away jet is not completely contained. 

If you refer to Appendix A of this note, you can convince about ESA allows rejecting events 
that some jet finders as JETAN could reconstruct as a perfect di-jet, although part of the jet 
stays out of the acceptance. 

We checked also the multiplicity distributions for the dijets events. In the Fig. 1131 we show 
the multiplicity distribution for the away and toward sides. 

A mono-jet event taken from region B of the ESA map is shown in the Fig. 1141 

The lines correspond to mctracks with pr > 0.8 GeV/c. For this event: 1 — T = 0.00216 
and R = 0.99874. Again, by counting primary charged particles within the azimuthal range: 
7r/4 rad < A(f> < 37r/4 rad, we estimate the total transverse momentum of the jet. In the Fig. 
1151 is the distribution of the total transverse momenta for events of region B. Note that the 
peak of this distribution is at ~ 6 GeV/c, in the case of di-jet events this peak is also located 
at ~ 6 GeV/c. 

The multiplicity distribution for mono-jet events is in Fig. 1161 

The conclusion of this part of the analysis is, that according with the results of the visual- 
ization and the behavior of the multiplicity and transverse momentum spectra ESA works fine 
for discriminating the di-jet events from the mono-jet ones. Using the event shape analysis 
the signals can be cleaned in order to improve the jet studies. 




Figure 12: Ratio of the transverse momenta of the toward jet over transverse momenta of the 
away jet as a function of R for dijet events. The participants have pt > 0.3 GeV/c. The near 
side corresponds to the azimuthal range: n /4 rad < A<j> < 37r/4 rad. For the away side: 5tt/A rad 
< A<f> < 77r/4. Charged component (right), (left) including the neutral component. 
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Figure 13: Multiplicity distribution of particles with transverse momenta: pt > 0.3 GeV/c for 
dijet events. The particles belonging to the leading jet (black line) consist of primary charged 
particles in the azimuthal range: 7r/4 rad < A<f> < 3ir/4 rad. The away side (red line) corresponds 
to particles with: 57r/4 rad < A</> < 77r/4 rad. The plot shows two cases: reconstruction (dotted 
line) and generation (solid line). 



Figure 14: Visualization ol a mono-jet event. The lines are mctracks with p t > 0.8 GeV/c. The 
arrows contain the particles associated with the jet. 




Figure 15: The distribution of the total charged transverse momenta of the identified mono-jets. 
The particles which were counted are within the azimuthal range: 7r/4 rad < A<f> < 37r/4 rad; and 
they have: p t > 0.3 GeV/c. The generated (solid line) and the reconstructed (dotted line) are 
shown. 
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Figure 16: Multiplicity of primary charged particles (p t > 0.3 GeV/c) of mono-jet events. The 
azimuthal range is: 7r/4 rad < A</> < 37r/4 rad. The results of generation (solid line) and the 
reconstruction (dotted line) are compared. 

4.2 Three-jet events 

The green distributions of the Fig. fallows to observe a double hump structure in the away 
side of the azimuthal distribution. They look like if the three particles with the highest pt in 
each event were distributed in the transverse plane according to: the leading particle at n/2 
radians, and the others at: 7n/6 radians and ll7r/6 radians, respectively. 

There are many configurations of isotropic distributions that could lead to high 1 — T 
values, however the scans of the events as shown in Fig. [17] for a typical event result in a 
clear three jets configurations. A small contribution from events with more than 3 jets can be 
found. 

It is important to say that this class of events occurs completely inside the acceptance of 
our detector, we have to remind that this is controlled by the term R. For example, if we 
increase the range of R, the azimuthal distribution in the away side shows a shift of the two 
peaks because the calculation of the variables use incomplete information of the event. 

In order to see the conservation of the transverse momentum in this class of events, we 
divided the azimuth into three regions: 

Near side: formed by particles in the azimuthal range: 7r/4 rad < A(f> < 37r/4 rad. 

Away side: formed by the particles in the remainder of the azimuth. 

As in the previous cases we have taken into account only primary charged particles with 
pt > 0.3 GeV/c. The Fig. [T5] shows the transverse momentum of each side. The agreement 
between both spectra is reasonable in the limits of the statistics. 

We compute also event by event the ratio of the transverse momenta of the toward jet over 
the vectorial sum of the away jets as a function of R. The result of this analysis is shown in 
the Fig. 1191 Again these distributions reach their peaks at 1, suggesting the correct transverse 
momentum conservation. 

The events with mercedes topology were found in MB simulation at 10 TeV in the c. m. 
as well as in MB simulations at 200 GeV |18| . 




Figure 17: Scan of one event of region C. Particles which appear in the picture have: pt > 0.5 
GeV/c. 
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Figure 18: Transverse momentum spectra for three-jet events, the including charged particles with 
Pt > 0.3 GeV/c. Toward side (black line) corresponds to particles in the azimuthal region: ir/4 
rad < A(j) < 3n/ 4 rad. Away side formed by particles in the remainder azimuth. The results at 
generator level (solid line) and reconstruction (dotted line) are shown. 
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Figure 19: Mercedes events. Ratio of transverse momentum distributions of toward jet over trans- 
verse momenta of the away jets. The A(f> ranges used are the same as in the previos figure. 



4.3 Multiplicity in the context of ESA. 

In this section we investigate the multiplicity of each event which we selected. In order to 
do this task we plotted the multiplicity spectra of the full sample (1 200 000 events), and we 
compared it with the multiplicities of the events with a given thrust value. The multiplicity 
is the number of primary charged particles in the acceptance \rj\ < 1 with p t > 0.3 GeV/c. 
As we can see in the plots of the Fig. 1201 the conditions which we demanded to each event 
reduces the number of low multiplicity events. The events with are not related with any of 
the classes we studied belong to multi-jet events, and also are of high multiplicity. We observe 
that the mercedes events are generally of large multiplicities. 



5 Bulk analysis 

The analysis of the event shape space projected on the 1 — T offers also interesting results. 
Without detailed analysis of tje jet components one may directly compare with the different 
generators. In the following we present the results of each analysis for the generators: Pythia 
and Phojet. Another interesting question is what happens with the events which thrust 
and recoil values are outside of the regions A, B and C. In the figure [21] there is the 1 — T 
normalized spectrum for different multiplicity bins. Three histograms are shown, the upper 
one corresponds to the cut R < 0.9, the middle: 0.9 > R > 0.4; and finally the events with 
R < 0.4 appears in the bottom one. It is clear that the distributions are quite different in 
the cases R > 0.9 and R < 0.4. In the last one the probability of finding a high multiplicity 
mercedes event is bigger than the probability of finding a low multiplicity event with mercedes 
topology. For events of region A, the maximum of the distribution is attained for the lowest 
multiplicity events. Note also, that in the middle R region (unexplored yet); the spectrums 
associated to different multiplicity bins are quite similar in their shape. 

5.1 Comparison between Pythia and Phojet in the context of 
ESA 

The Monte Carlo generators Phojet [19] and Pythia [20] use both LO QCD matrix elements 
for the hard scattering sub-processes. Initial and final state parton radiation and the string 
fragmentation model are included as implemented in the JETSET program [21] , The two 
Monte Carlo generators differ in the treatment of multiple interactions and the transition 
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Figure 20: Multiplicity (charged paricles with p t > 0.3 GeV/c) distribution. In black is the full 
true spectrum. The different colored distributions correspond to the three classes of topologies 
which we discriminated using ESA. 



from hard to soft processes at low transverse parton momentum. The hard parton-parton 
cross-section diverges towards low transverse momenta and therefore needs a regularization to 
normalize to the measured total cross-section. 

Hadronic collisions at high energies involve the production of particles with low trans- 
verse momenta, the so-called soft multi-particle production. The theoretical tools available at 
present are not sufficient to understand this feature from QCD alone and phenomenological 
models are typically applied in addition to perturbative QCD. The Dual Parton Model (DPM) 
[22] is such a phenomenological model and its fundamental ideas are presently the basis of 
many of the Monte Carlo implementations of soft interactions. 

The Monte Carlo event generator Phojet can be used to simulate hadronic multi-particle 
production at high energies for hadron-hadron, photon-hadron, and photon-photon interac- 
tions with energies greater than 5 GeV. It implements the DPM as a two-component model 
using Reggeon theory for soft and leading order pQCD for hard interactions. Each Phojet 
collision includes multiple hard and multiple soft pomeron exchanges, as well as initial and 
final state radiation. In Phojet pQCD interactions are referred to as hard Pomeron exchange. 
In addition to the model features as described in detail in [T5], the version 1 .12 incorporates 
a model for high-mass diffraction dissociation including multiple jet production and recursive 
insertions of enhanced pomeron graphs (triple-, loop- and double-pomeron graphs). 

So, Phojet provides an alternative to Pythia for the study of processes that cannot be 
calculated with pQCD, such as minimum bias events (events with high cross section and low 
transverse momentum) and the underlying event activity in events with a high transverse 
momentum parton-parton collision. 

5.2 Implementation of ESA 

The Pythia and Phojet events used are part of the standard simulations in ALICE. The 
samples are staged on alicecaf and they correspond to proton-proton collisions at 10 TeV in 
the c. m., minimum bias, magnetic field of 0.5 T. This analysis was done using 200 000 events 
in each sample. 

The Fig. !22l shows the comparison between the normalized spectrum 1 — T for the samples 
generated with Pythia and Phojet. Three different regions in R are explored: the first one: 
populated by single-jet events (R > 0.9). The second is the intermediate R zone (0.4 < R < 
0.9). And finally the zone where there are di-jet and three-jet events (R < 0.4). 

Note that in agrement with the discussion of the previous section, we can observe in the 
region of the highest values of R and 1 — T a more copious presence of mono-jet events in the 
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Figure 21: Projection on the 1-T axis as a function of R and multiplicity. The regions with 
interesting topologies are indicated. 
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Figure 22: Projection on the 1-T axis for Pythia and Phojct generated events. 



simulations generated by Phojet than in the case of Pythia ones. The same situation appear 
in the case of di-jet events. Also, note that the particles generated by Pythia look distributed 
in a way more isotropic compared to the Phojet sample. For the more isotropic topologies and 
especially the mercedes events the production with Pythia is more copious than with Phojet. 



6 Conclusions 

In the present note we have demonstrated the applicability of the Event Structure Analysis 
in the case of measurements of charged particle tracks in a detector of limited acceptance like 
ALICE. The phase space in recoil and thrust variables offers a wide variety of possible uses: 

Selection of well characterized events like monojets, dijets and three jets ones. ESA is 
in that instance not a replacement for any kind of jet finder but more of a preselection 
that can be then studied with jet algorithms. As we report the number of "clean" dijets 
for instance is minute in comparison with the total sample. The use of conventional jet 
finding algorithm in events with many jets as is the case at LHC leads inexorably to the 
use of cuts that are sometimes leading to results difficult to interpret. 

Bulk comparison of the event characteristics with existing models. Since in the majority 
of cases we are confronted with multijets events we believe useful to establish a way to 
compare model predictions with projection on the sphericity axis of events belonging to 
different recoil and/or multiplicity intervals. 

Last but not least any deviation of the present predictions of generator would lead to a 
rapid detection of events located in an "odd" region of the R vs. (1 — T) phase space. 

We have presented the discrimination power of ESA for specific topologies: dijets events, 
monojets and three jet topologies. According to the results presented, the ESA may be of 
use in the physics analysis starting with the first data since with a sample of 400000 events of 
minimum bias events, there are many events with the topologies which we discussed. 
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Figure 23: Analysis of the slice 1 — T < 0.05 of the thrust map. The jet finder UA1 was used for 
reconstructing the jets. 



A Appendix 

A.l About the usefulness of ESA in combination with JETAN 

The number of jet to be encountered in pp collisions at the LHC will be very large. In the 
general case we will find many overlapping jet cones as is shown in Fig[T] On the other hand 
the small acceptance for jets causes that many jets are only partially in the acceptance. We 
illustrate that with Fig |23l where we present the ratio of the difference of the total transverse 
momenta found for events identified by JETAIS0 as dijets but belonging in the ESA analysis 
to two different slices in R. For the lowest one R < 0.2 we see that the difference of momenta 
is generally small while for the slice of the highest R the difference is markedly broader. We 
therefore believe that ESA apart from other virtues for the analysis of bulk properties of 
the events has a very important task to play in the identification of specific topologies. The 
subsequent comparison of these identified topologies with any kind of jet finder allows in our 
mind to study "well cleaned" samples without resorting to different cuts as is usual in the use 
of the jet finders. To illustrate our point of view we show two figures: Fis l24l the distribution of 
particles in the eta phi acceptance of ALICE using events identified as dijets by JETAN, and 
on the right side the distribution obtained for events identified as dijets by JETAN applied to 
ESA dijet events. A completely clean separation between the near and far side jets is visible. 
However, applying JETAN to the ESA dijet sample we found that a number of events were 
identified in the present use of ESA as monojet events. The reason is clear from the eta phi 
distribution in Fig |25l where we plot the "monojet" events found by JETAN. The reason is 
simple - in ESA we did not limit the eta range for the leading particle of the "away" jet. Hence 
a part of the events are de facto encountered in the edges of the eta acceptance. The detailed 
study of these distribution will be pursued. 



6 JETAN is a module of AliRoot which includes different jet finders. In this analysis we used the UA1 jet finder 
algorithm which is based on a cone type [23] ■ The cut in the transverse momentum of the particles which we included 
is: pt > 1.5 in the cone radius: %J (A</>) 2 + (Arj) 2 = 0.7. And the minimum transverse energy of the jets: 1.5 GeV. 



Figure 24: Topologies of the events identified as dijets. The transverse momentum distributions 
of the primary tracks in the plane rj vs A<j) (the leading particle is at Aip = tt/2) for different 
samples of events are shown for the following cases. Dijets reconstructed with JETAN (left). 
Dijets identified by ESA and JETAN (right). 



2J ESA, 1J JETAN 




Figure 25: Topologies of the events identified as dijets by ESA but reconstructed as monojet by 
JETAN. Each bin in rj-A(cf> — 7r/2) is weighted with the transverse momentum of the charged 
particles. 
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